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Abstract 

We consider the spin 1/2 model coupled to a slowly varying magnetic field in the 
presence of a weak damping represented by a Lindblad-form operators. We show 
that Berry's geometrical phase remains unaltered by the two dissipation mechanism 
considered. Dissipation effects are twofold: a shrinking in the modulus of the Bloch's 
vector, which characterizes coherence loss and a time dependent (dissipation related) 
precession angle. We show that the line broadening of the Fourier transformation 
of the components of magnetization is only due to the presence of dissipation. 
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1 Introduction 



The existence of geometric phases in quantum systems, ever since its discovery by Berry 
H|, has attracted considerable interest both from the theoretical and experimental view- 
points. Several applications of this phenomenon in different areas of physics have also 
been studied^, [| |J. Comparatively much less work has been devoted to the question of 
the dynamical evolution of such systems in the presence of a weakly dissipative medium. 
Basically the available results can be summarized as follows: nonhermitian operators lead 
to a modification of Berry's phase^], §, [7|], stochastically evolving magnetic fields produce 
both energy shift and broadening ||, phenomenological weakly dissipative Liouvillians 
alter Berry's phase by introducing an imaginary correction^] or causing damping and 
mixing of the density matrix elements |IIJ. Ellinas et al. || obtain their results study- 
ing the eigenmatrices of the complete Liouvillian superoperator in a time-independent 
basis, while Gamliel and Freed |10[, find closed formal expressions for the density matrix 
representation in the instantaneous frame of the hamiltonian in the adiabatic limit, and 
in the weak dissipation approximation. However, even in this regime the results in the 
cases studied can only be extracted through approximations (in addition to the adiabatic 
approximation) or numerical computation. 

In the present work, we consider the celebrated example of a spin 1/2 particle coupled 
to a slowly evolving magnetic field in the presence of a weakly dissipative medium as 
represented by a Lindblad form superoperator JTT], [12], [13|], incorporated in a convenient 
and physically motivated frame. In the absence of dissipation, a simple geometrical in- 
terpretation of the results emerges in terms of both the Bloch's vector || and the phase 
vector JT4j. The geometrical phase appears as a "delay" or an "advance" in the precession 
period of the Bloch's vector with respect to the period dictated by the magnetic field's 
frequency. The precession of this vector occurs at a fixed angle with respect to a fixed axis 
about which the external magnetic field precesses, as usual. We will introduce the dissi- 
pation via a semi-group type dynamics, such that dissipation does not alter the precession 
frequency. However, the x, y and z components of Bloch's vector are altered in different 
ways by nonunitary effects. In particular, the modulus of Bloch's vector will shrink. The 
relation between the reduction of Bloch's vector modulus and the loss of coherence has 
been explored by Stodolsky and collaborators |T^, |H| [TTJ. In order to make it explicit, we 
introduce the linear entropy (or idempotency defect) as a measure of purity lossfl8[ 
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S(t) = l-Trp 2 (t) (1) 



and write 




S x (t) - iS y (t) 
1 - S z (t) 



= -(l + S(t)-a), (2) 



6(t) 



l-\S(t)\ 2 



(3) 
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which confirms that the shrinking of the modulus of S(t) is a measure of coherence loss. 

In order to discuss dissipative effects on interference patterns due to Berry's geo- 
metrical phase the natural framework is to consider the evolution of the density matrix, 
which completely characterizes the interference effect due to the different geometric phases 
acquired by the eigenvectors of the hamiltonian of the system, even in the absence of 
dissipation 0. Of course, this time we add a nonunitary Liouville operator contribution 
to the dynamics. 

In section 2 we present the two level quantum system of a spin 1/2 in the presence of 
an external magnetic field precessing with constant angular velocity around a fixed axis. 
The master equation of the model is written in what we call the diagonal frame where 
Lindblad superoperators are introduced to describe the nonunitary part of the Liouville 
operator. In subsection 2.1 we consider the adiabatic limit of the quantum system in 
thermal equilibrium with a reservoir of electromagnetic fields and in subsection 2.2 we 
consider the case of a dephasing process. In both cases, we treat the coupling of the 
quantum system with its environment in the weak regime. In section 3 we show how the 
geometric phases and the dissipation effects acts on the z component of the magnetization 
of the spin 1/2 coupled to a reservoir of electromagnetic fields at thermal equilibrium. In 
section 4 we summarize our conclusions; finally, in appendix A, we present the non-unitary 
part of the liouvillians of the models under consideration in the instantaneous basis of the 
hamiltonian. 
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2 The adiabatic limit of the spin 1/2 model in the 
weak coupling regime 

We consider a spin 1/2 variable (two level model) coupled to a time dependent magnetic 
field precessing around the z-axis. The unitary contribution for this evolution is given by 
the hamiltonian 

H s (t) = na-B{t) 

/cos(£) sin(0)e--*\ 
^ V sm^je*"* - cos(#) J ' 1 > 

written in the basis of the eigenstates of the z-component of the spin, where B is the 
norm of the external magnetic field, 9 its azimuthal angle, uo the precession frequency and 
the constant p = — being g Lande's factor and p R the Bohr magneton. We are using 
natural units (c = h = 1). For the sake of later calculations, it is convenient to define 
two unitary transformations: the first one, R(u>,t), takes us to the rotating frame where 
the hamiltonian is no longer time dependent; the second one, D(B,8,u), diagonalizes 
the effective hamiltonian (time independent) that drives the dynamics of the final matrix 
representation of the density operator. After the first transformation, R(u,t) = 
the density matrix and the hamiltonian read 



-<7z 



p R (t)=e^p(t)e-^% (5) 
and 

H R = fiB{sm{6)a x + cos(6)a z ). (6) 

In analogous manner, after the second transformation we get, in the diagonal frame, the 
density matrix 

p D (t)=B T p R (t)B, (7) 

and the effective hamiltonian 



H D = B T (H R - |<7 2 )D 
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where Ai = \ p 2 B 2 sm 2 (9) + (pB cos(9) — tt) 2 . The rotation matrix D = D T is equal to 



D = \l\ - ^-(/xB cos(e) - |) a x + ^ + ^-(/*Bcos(0) - |) (9) 



One possible way to add dissipative contributions to the above dynamics is to include 
a Lindblad type superoperator in the evolution equation in the diagonal frame 



j t PD (*) = -z [Aia 2 , p D (t)] + kC D p D (t), (10) 

where is the dissipation constant. The weak coupling regime is characterized by the 
condition 1. 

Our aim is isolating, in the density matrix, the effects of dissipation on interference 
due to geometric phases. The representation of the density matrix in the diagonal frame 
is not very enlightening for this purpose. This is better realized in some basis of the 



instantaneous eigenvectors of hamiltonian @, as was done by Gamliel and Freed |10| . We 
define pi(t), the matrix density in a basis of the instantaneous eigenvectors of hamiltonian 
(ID EH- The relation between pi(t) and poif) is 



Pl (t)=V\t)Vp D (t)VV(t), (11) 
where the matrix V(t) is equal to 



vw;=, . jK 2 " Sm( !lfi I • (12) 



sm^|je 2 cos^je 2 



The time evolution of P i(t) is given by 
d 



where 



— pj (t) = -i (jiB + t£)<Tz- 7jVn(t),pi (t) + kC T pi(t) (13a) 



a(t)-[ C ° SW -M0)e~^ . , . 

an[t) ~ I -sm(0)e^ -cos(0) > 1 ' 
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and Cipi{t) is obtained from CoPoit) through a similarity transformation equivalent to 

Before we specialize our discussion to any particular liouvillian, we study the adiabatic 
limit of eqs. fll3|) for the coupling constant k in the weak regime. In this regime, the matrix 
Cipi(t) is written as a linear superposition of the elements pjj(t). We remind that the 
density matrix of a two level model must satisfy two conditions: i) Tr(pj(t)) = 1 and 
= (Pn (*))*■ As a consequence of those conditions, the density matrix has only 
two independent elements. We take the elements p{i(t) and p[ 2 {t) as our two independent 
entries. The general form for the time equations of those two elements in the weak 
coupling regime is 



^Pn(t) = a n (uj, k; t)p[ x {t) + a 12 (co, k; t)p[ 2 {t) + a n (u, k; + b x {u, k; t) (14a) 



and 



^pLW = fl 2i(^, k; t)pn(t) + {-2ipB + a 22 (u, k; t))p[ 2 (t) + a 23 (uj, k; t)p 21 (t) + 

+b 2 (u,k;t). (14b) 

We make the change of variables: 



Ptjit) = e-^-^ p tJ (t) (15a) 

with Ei = fiB and E 2 = —p,B. We may disclose the time scale T in the differential 
equations by introducing the following transformation upon the time parameter P0|, PT|: 



s = |, (15b) 



where T — — . 

With the new variables, eqs.(|T4]) become 



^-Pix{s) = Taii(u,k;s)pii(s)+Tai 2 (u,k;s)e- 2l » BTs pi 2 (s) + 
as 

+ Tai 3 (co,k;s)e 2iflBTs p 2 i(s)+Tbi(to,k;s), (16a) 

^-Pi 2 (s) = Ta 2 i(u, k; s)e 2 ^ BTs pn(s) + Ta 22 (u } k; s)p 12 (s) + 
as 

+ Ta 23 (co, k- s)e^ BTs p 2 i{s) + Tb 2 (co, k; s)e 2 ^ BTs . (16b) 
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We take the differential equation for pi 2 (t) to exemplify the discussion of the adiabatic 
limit of eqs.([16|). At this point we will follow closely the references [[H], pi]] . We point out 
that the adiabatic approximation is not recovered by an uj expansion of the terms on the 
r.h.s. of eqs.(|I"6l). 

The Volterra equation|22| obtained from eq. (|16b| ) is 



Pl2{S) 



ps ps 

= p 12 (0)+T/ a 22 (cu,k;s')p 12 (s')ds' + T a 21 (co, k; s')e 2i » BTs ' p n {s')ds' + 
Jo Jo 

+ T [ S a 23 (cu,k;s')e A ^ BTs ' p 21 (s')ds' + T f b 2 (cu, k; s')e 2 ^ BTs ' ds' . (17) 
Jo Jo 



In the limit T — > oo, the Riemann-Lebesgue Theorem [£3[] gives that 



T->oo 



lim / F(s')e tailTBs 'ds' = 0, 



if F(s') is a piece-wise continuous function in the interval [0, s] and a G R. As a conse- 
quence of this theorem the last three terms on the r.h.s. of eq.(|l7]) vanish in the adiabatic 



limit ->• T -»• oo). 



Integrating by parts eq.(|18|) we get 



/•(•s')r "'"' /7 -S/,' = i-L F(s)< : """- - F(0))] - 4t^7T / ^7 (/••(>')) e^^ds'. 



iT apB 



iT apB 



(19) 



From eqs.([L3|) and ( p.4b|) , the coefficients a 2j - and b 2 have the general dependence on u 
and k: 



a 21 (u,k;s) = U a^(s) + k[a^(s) + ^A 21 (s)+oU^] 



a 22 (uo,k;s) = - cos(6)) + k [ af 2 \s) + -^A 22 (s) + O 

a 23 (u,k;s) = k[d { S(s) + ^A 23 (s) + oU-^j 



oo 



and 



oo 



\ 



pB J 



b 2 U'. k: .s) = oob^(s) + k [ bf\s) + ^B 2 (s) + O I I 



fiB J 



(20a) 
(20b) 
(20c) 

(20d) 
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The integrands of all integrals on the r.h.s. of eq.flTTD which contain the oscillatory 
function e ia ^ BTs ' ( with a = 2 or 4, after integration by parts, each of those integrals has 
its order in T decreased by one unit, and acquires a multiplying constant of value 
or jjj. From the conditions satisfied by the two simultaneous regimes: i) adiabatic limit 
(-^ < 1) and ii) weak coupling limit (-^ < 1), we can neglect those terms in comparison 
to the first two terms on the r.h.s. of eq.([[7]). The two previous inequalities do not impose 
any constraint to the ratio though. The differential equation satisfied by pn(s) in the 
adiabatic limit and weak coupling regime is 



d 
ds 



Pl2{S) 



T 



— cos(#)) + ka^is) 



Pl2{S) 



(21) 



The term proportional to ^22(5) was dropped, since it is of higher order in (j^J- For 
k ^> ou, the off-diagonal elements of the density matrix vanish before the external magnetic 
field B(t) returns to its configuration at t = 0. The other uninteresting situation from the 
point of view of the appearance of an imaginary correction to the geometric phase is the 
condition k <C uj when the dissipation effects can still be neglected after one period. In 
this work we discuss the case k ~ u when both terms on the r.h.s. of eq.(^Tj) contribute 
to the time evolution of pi2(s). 

By a similar discussion we obtain the time equation of pu{s) in the adiabatic and 
weak coupling regimes, 



d_ 

ds 



pn(s) = Tka ( S(s)p 11 (s) + T ub^>(s) + kb { ?\s 



■Mi 



(22a) 



From eqs.fllq) we can affirm that 



a n (uj,k;s) = k[a^(s) + ^A n (s) + oU^J )), (22b) 



h{uj,k;s) = k[bf\s) + ^-B 1 (s) + ol (^-) 



pB 



\pB J 




(22c) 



In order to verify if we can get imaginary phases from eqs. QzTp and (|22a[ ) due to the 
coupling of the quantum system to a dissipative medium, we write the density operator 
at t = in the instantaneous basis of hamiltonian (f|) 
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2 

P(0) = E «ii(0)|#W?|. (23) 

We have H(O)|0j) = #j(O)|0°), j = 1 and 2. Being the time evolution of the 

eigenvector we have 

2 

p(*)=x; ^(t)i^°(t))(^(t)i, (24) 

3,1=1 

where the time evolution of is driven by H(t). Differently from Gamliel and Freed, 

we include the phases coming from the unitary evolution and the geometric phase in the 
dyadic product \4>^(t))((p^(t)\. In the adiabatic approximation, we get 

\<fjj(t)) = e^' w e-^^lfat), (25) 

where jj(t) is the geometric phase acquired by the eigenvector \4>j',t) is the instan- 
taneous eigenvector of H(t) ( H(t)| = E^t^ft) ) and = f /„* dt'E^t'). 
The density operator at any time, in the adiabatic limit is 

2 

p{t) = <*ji(t) j^V-vW e -*«^(*)>-<^iC*)»*|0 i; t|. (26) 
j,i=i 

From eq. (|26|) we recognize that the phase e l in the element p^if) is just the 

difference of the geometric phases of the instantaneous eigenstates \<pj] t) and \<f>i\ t) in the 
absence of dissipation. 

The dynamics of the coefficients atji(t) is ruled by the nonunitary evolution of the quan- 
tum system and it is independent of the particular choice for the instantaneous eigenstates 
of the hamiltonian, up to a multiplicative constant^]. In the model under consideration 
(see hamiltonian (Q)), the eigenvalues Ej(t), j '• = 1 and 2, are time-independent. It is 
simple to get the time equations of an(t) and au(t) from eqs. ( |22a| ) and ([H]), respectively 

fc(aSJ(t)a u (t)+6f )) (t)) J (27a) 
ka { S[t)a l2 {t). (27b) 
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d . . 

d . . 
dt ai2{t) 



The constants a^(t), o^(t) and bf\t) depend on the particular master equation that 
describes the behaviour of the quantum system interacting with the dissipative medium. 
In the next sub-sections we consider two particular interactions of the two level model with 
a reservoir: i) two level model in thermal equilibrium with a reservoir of electromagnetic 
fields; ii) dephasing process in a two level model. 

2.1 Adiabatic limit of a two level model in thermal equilibrium 

As discussed before, to incorporate the dissipative effects in the two level model we intro- 
duce the Lindblad superoperator in the diagonal frame. The master equation of the spin 
1/2 model coupled to a reservoir of electromagnetic fields in thermal equilibrium in the 



diagonal frame is [111], E5 



-j t PD (t) = -i [Aicr 2 , p D (£)] +k(n+l) {2a_p D (t) a + - p D (t) a + a_ - a + a_p D (t)) + 
+kn (2a + p D (t) cr_ - p D (t) a_o + - o_o + p D (t)) , (28) 

where k is the dissipation constant at zero temperature and n is the average number 
of excitations of the weakly coupled thermal oscillators at inverse temperature (3. An 
important requirement for the introduction of this Lindblad type superoperator is that it 
leads asymptotically to a thermal equilibrium. In appendix A we give the master equation 
of this physical process in the instantaneous basis of hamiltonian (|4]) for arbitrary value 
of to. 

From the master equation in the instantaneous basis of the hamiltonian we obtain the 
equations for «n(t) and ai 2 (t). These equations in the adiabatic and in the weak coupling 
limits become 



_a n (t) = -2k(l + 2n)an(t) + 2kn, (29a) 

j t a 12 (t) = -k(l + 2n)a 12 (t). (29b) 

The constant k does not come up on the r.h.s. of eqs.(|29|) due to the time variation of any 
classical parameter that characterizes the reservoir. From its explicit definition [[11], [25{] 
it can not be written as: f(t)g(t), where f(t) and g(t) are two regular time dependent 
functions. The same is true for n. Therefore the imaginary phase x{t) defined as: au{t) = 
a 12 (0)e ix ^\ with 
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x(t) 



t 

i, I k(l + 2n)dt' 
'o 



(30) 



is not geometric. 



The solution of eq. (|29a|) is 



a n (t) 



n 



l + 2n 



+ 



<*ii(0) 



n 



2n 



-2k f*(l+2n)dt' 



(31) 



The exponential decay on the r.h.s. of eq.([H|) means that the population of the 
instantaneous eigenstates of hamiltonian (|J) varies in time and consequently the Adiabatic 
Theorem is not valid for this dissipation mechanism. 

Exactly soluble models are always important checks to approximation schemes. Eq . (^8|) 
is exactly solved and the solutions are 



n 



2n + 1 



3 -2fc(2n+l)t 



+ Pn(0)e- 



-2k{2n+l)t 



(32a) 



(32b) 



It is straightforward to obtain the adiabatic and weak coupling limit eqs. ( |29l) from eqs. 032| 



2.2 Dephasing process in two level system 



Another interesting process well studied in the standard textbooks [|11], [25|] is the phase 
destroying process which might appear due to elastic collisions. In general, those effects 
are incorporated in the master equation of the two level model, besides the energy dissi- 
pation process studied in subsection 2.1. Since we are studying the coupling of the spin 
1/2 to a dissipative medium in the weak coupling limit, the inclusion of the phase de- 



stroying process in eq.(|28|) gives corrections to the coefficients dij(uj } k; t) in eqs. (P9|) . Due 
to the linearity of the equations, the new imaginary phases coming from the dephasing 
process are added to the ones obtained previously. For the sake of simplicity we study the 
imaginary phases acquired by the variables an(t) and a^it) only due to the dephasing 
process. The master equation written in the diagonal frame is 
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d k 

—p D (t) = -i [\ia z , p D (*)] + - {o z p D {t)o z - p D {t)) . (33) 

In appendix A we present the master equation of the spin 1/2 with the dephasing 
effect included for arbitrary value of the angular velocity to of the external magnetic field. 
In this subsection we study the time equations of coefficients au(t) and ai 2 (t) in the 
adiabatic and weak coupling regimes. Taking into account our discussion in section 2, 
eqs.( |A.2[ ) and ( |2"7| ) we obtain 



^an(t) = 0, (34a) 
at 

j t a 12 {t) = -ka 12 (t). (34b) 



that have the solutions: 



an (t) = an(0), (35a) 
a 12 (t) = a 12 (0)e-ti kdt '. (35b) 

From eq. ( |35a| ) we conclude that the Adiabatic Theorem is valid in this process, since the 
population at each quantum state does not vary along the adiabatic process. By analogous 



reasons to the ones discussed in subsection 2.1, the phase in eq . (|35b| ) is not geometric but 
a time dependent imaginary phase that destroys the off-diagonal elements of the density 
matrix. Eq.(^) is exactly solved. The solutions ( pop are easily recovered from the exact 
solutions when we calculate them in the adiabatic and weak coupling limits. 

In order to understand why references f|-[[7|] give an imaginary correction to the ge- 
ometric phase and the models studied here do not, we compare eq.(2.3) of reference || 
with our eqs. (|3~0D and (|35b| ). The coefficient that multiplies of variable C\(t) in eq.(2.3) 
of reference has the form (0{t)\4:\if)(t)). Since the vector states depend on a periodic 
external parameter R(t), this coefficient, in the adiabatic approximation, corresponds to 
a correction to Berry's phase written as a closed curve in parameter space. The coeffi- 
cient k(l + 2n) that multiplies a\ 2 (t) in eq.( |29b| ) does not arise from any variation of an 
external parameter. The same is true for the coefficient k in eq. ([34b| ). That is the reason 
that allows us to claim that the exponentials acquired from eqs.(^) and (|35b|) have no 



geometric origin, which means that the suppression terms are a function of time instead 
of some path parameter. 
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3 Contribution of the geometric phase to the 
magnetization 



In order to verify how the geometric phases and the dissipation affect the physical quan- 
tities, we return to the spin 1/2 model in the presence of a reservoir of electromag- 
netic fields in thermal equilibrium. Under the initial condition p(0) = |^(0))(^(0)|, with 
\ip(0)) = cos(a)|+) + sin(a)|— ), in the adiabatic limit and weak coupling limit we end up 
with 



pf{t) = i sin(2a - 9) e" 2 ^* e ~ k{m+1)t e -<«-(i-««(»))*. (36b) 

We have introduced the Berry's phase tracer a that helps us identify the contribution of 
Berry's phase to physical quantities. At the end, we take the tracer equal to one. 

Thus we see that, in this model, dissipation does not affect Berry's geometrical phase, 
but only makes it harder to observe their interference effect: such information is contained 
in the off-diagonal terms of the density matrix in the instantaneous frame, which vanish. 
By the way, we can observe that the system considered in this contribution is analogous to 
the case of the classical Foucault pendulum where dissipation diminishes the amplitude, 
but do not affect the rotation of the oscillation plane. 

The definition of the geometrical phase can also be given in terms of a phase vector as 



in reference ||14j|, where the discussion is confined to pure states. It is, however, a relatively 
simple matter to extend the definition of the phase vector for mixed states. In this case 
we find that the loss of coherence will shorten the phase vector in a manner which is 
completely analogous to what happens to Bloch's vector. This norm reduction, as in the 



present case, reflects the asymptotic vanishing of off-diagonal density matrix elements p6 

Since we are interested on interference effects due to geometric phases, Bloch's vector, 
defined in eq.(^), is suitable to provide a graphic visualization of the density matrix ||. 
From eq.(0), we have that these effects are clearly described by Bloch's vector in the 
instantaneous frame. In this frame the projection of Bloch's vector sweeps the xy plane 
and makes an angle smaller than that of the magnetic field by an amount proportional to 
the solid angle Q(9), while its length decreases exponentially with a time rate of k{2n + 1). 
We illustrate in figure [l] the shrinking of the projection of Bloch's vector in the xy plane, 
due to the presence of dissipation. Even though we have discussed the case k ~ u, 
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we choose k/u = 10 in order to show more clearly the decreasing of this projection 
in an interval An/fiB. Shortening in z is faster than in the xy plane, causing a time 
dependent azimuthal angle. Figure [2] shows the plot of the time evolution of Bloch's 
vector for k/u = 0.2. For t —>■ oo Bloch's vector has only non-zero ^-component and 
S z {t -> oo) = . 

For time t, such that t ^ nT (T = 2ir/u and n is an integer) the off-diagonal elements 
of the density matrix pi(t) depends on the chosen condition satisfied by {4>i\ t\^\<j>%\ t) fT9fl . 
For the sake of comparison with experiments it is necessary to calculate dissipation and 
geometrical phases effects on measurable quantities. In this model the natural candidates 
are the components of the magnetization vector (m)(t). Let us consider the ^-component 
of magnetization whose average value has the expression Tr(p(t)m z ). Since the trace is 
independent of the particular basis applied to calculate it the result is independent of our 
particular choice of t). In the adiabatic approximation and the weak coupling 

limit, using eqs.fl36|), we get the Fourier transform of (m z )(t), 



(m z )(u') = -£= 7t8{u') - i cos(#) ( 7 ^- T + cos(# - 2a)) ; ; - L_ . _ + 



L2n + 1 v ' w V2n+l v ') u' + 2ki{2n + 1) 

—i - sin(60 sinf^ — 2a) ( — H — -) , 

2 v; v ' \uj' + Y + ki{2n + l) w' - T + ki(2n + 1) / J' 

(37a) 

where a is given by the initial condition, and the resonant frequency T is equal to 



r = 2fiB - aucos(9). (37b) 

The first term of eq. (|37a|) corresponds to the constant component of the magnetic field, 
and the second one shows the dissipation effects on this field component. The last term 
on the r.h.s. of the above equation displays a real frequency shift which contains the 
contribution of the geometrical phase, as can be seen from eq. (|37b| ), due the presence of 
the tracer a, and a line broadening caused only by the dissipative evolution. These effects 
on the magnetization agree with the ones derived in reference where the path integral 
formalism was applied. 

The expressions for the other components of magnetization are analogous, but some- 
what lengthy. 
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4 Conclusions 



In summary, we have presented an analytical solution of the adiabatic limit of a spin 1/2 
in a precessing magnetic field embedded in a weakly dissipative medium, introduced phe- 
nomenologically. We consider two distinct nonunitary contributions that were accounted 
for by a Lindblad type superoperator in the diagonal frame. We are able to derive ana- 
lytical expressions for the geometric and imaginary phases in both cases in the presence 
of a weak dissipation in the adiabatic limit without further approximations. 

From eqs.(|30D and ( |35b|) we get that the nature (path dependent or time dependent) 
of the imaginary phase acquired by ot\2{t) depends on the mechanism that introduces the 
dissipation in the quantum system. In both cases that we have studied, the dissipation 
is present due to the two level system being in contact with a reservoir. The constant 
k on the r.h.s. of eq. (|30]) is associated to the time rate of population and not due to 
the variation of any external parameter. Consequently the quantum geometric phase for 
k = is not modified by a complex value. An analogous argument is valid to explain 
why the imaginary phase acquired by the entry p[ 2 (t) in the dephasing process is not 
geometric, either. In this last model the Adiabatic Theorem continues to be true while 
it is not true anymore for the spin 1/2 coupled to the electromagnetic fields at thermal 
equilibrium. Differently from Ellinas et al. in reference || we do not call this complex 
phase as Berry's phase. We reserve the name of "Berry's phase" only to phases (real or 
imaginary) that are path dependent. 

Decoherence effects are present and their manifestation is the shortening of the three 
components of the Bloch's vector. The fact that the dissipation effect causes the sup- 
pression of interference patterns due to the geometric phase is not a particular result for 
the chosen liouvillian. It is rather general, stemming from the fact that the dissipation 
mechanism is not related to the variation of any set of external periodic parameters. 
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Appendix A: Master equations in the instantaneous 
basis of the hamiltonian 



The master equation fl28D in the instantaneous basis of the hamiltonian for arbitrary value 
of oj is 



-nPi (t) = -i {v B + |) °z - ^n(t),Pi (t) k£, P , ( / ) 



(A. la) 



where a n (t) is given by eq. (|13b|) and 



C lPl (t) 



2n + l 



- 2 Pl (t) + (1 - A ) a n (t) Pl (t)a n (t) - e~ 2lM a + {t) Pl (t)a + {t) - 



e 2tuJt a.(t) Pl (t)a4t) + (1 + A ) a + (t) Pl (t)a^t) + a„{t) Pl {t)a + {t) 



- AVT^[e iuit (a n (t) P i(t)a-(t) + a_(t) Pl (t)a n (t) 



+ e- iujt (^ + (t) Pl (t)a n {t) + a n (t) Pl (t)a + (t)^ 



+ 



-{{ Pj (jt), Aa n (t) + VT^(e- luJt a + (t) + e^V_(*)} + 2A a + {t) Pl {t)a-{t) 



a. 



.(t) Pl (t)a+(t) 



e-^[a n (t) Pl (t)a + (t) - a + (t) Pl (t)a n (t) + 



+ e^(a_(t)p/(tK(t) - ff n (t) P/ (t)cr_(t))] }. 
In eq.( |AT5| ), we define: A = ±(^ P Bcos(9) - f J and 



(A.lb) 



a + (t) 



Suit 



sin ( g ) „ J _f nna r e \\2 e -i U t„ 



(sin(^)) 2 e-V. 



(A.lc) 



and er-(t) = ^cr + (t) 



The master equation (j33|) in the instantaneous basis of the hamiltonian for arbitrary 
value of oj is 
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d 

Jt' 

where 



PI (t) = -i [(jiB + |) a z - ^a n (t), Pl (*)] + \CiPi{t) (A.2a) 



£jp,(t) = A 2 a n (t)p I (t)a n (t) + AVl - A 2 [e^* (<7 n (t)p/(*)<7+(*) + a + (t)p 7 (t)a„(t)) + 
+ e^fa^p^V- (*) + <7-(*)P/(*K(*))] + (1 - A 2 ) [e- 2i ^a + (t) Pl (t)a + (t) + 
+ e 2 ^V_(t)p 7 (t) ( r_(t) + <T + (t) Pl (t)<T-(t) + <r_(*)pj(*)<T+(t)] - Pi(t). (A.2b) 
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Figure 1: Plot of the projection of Bloch's vector in the xy plane. We take a = 0, 6 = j, 
uj I jiB = 1CT 3 and k/ jiB = 10 -2 . Time evolution is plotted in the interval t = [0, -^]. 



Figure 2: Curve described by Bloch's vector in space. We take: a — 0, 9 — |, tu/fiB = 
1CT 3 and k/fiB = 2.10" 4 . The vector starts with positive z-component but ends up on 
the negative z-axis, after a time interval equal to the period of the external magnetic field. 
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